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Abstract. We introduce a method to investigate the longitudinal momentum spread 
resulting from strong-field tunnel ionization of Helium which, unlike other methods, 
is valid for all ellipticities of laser pulse. Semiclassical models consisting of tunnel 
ionization followed by classical propagation in the combined ion and laser field 
reproduce the experimental results if an initial longitudinal spread at the tunnel exit is 
included. The values for this spread are found to be of the order of twice the transverse 
momentum spread. 
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1. Introduction 

Applying a strong electromagnetic field bends the binding potential of an atom or a 
molecule, allowing the electron to tunnel out. In attosecond science PQ, this tunnel 
ionization is commonly assumed to be the starting point of many important phenomena, 
including High Harmonic Generation (HHG). Consistent semiclassical models are 
required to interpret and explain results from ultrafast laser experiments [2j [3] as well 
as help design new experiments. 

An assumption that comes from the tunnelling limit of the Strong Field 
Approximation (SEA), is that at the tunnel exit, an electron does not have any 
momentum parallel to the electric field (see for example [4]). This is contrary to the 
transverse momentum spread, for which the ADK theory [SJ El Ej predicts a Gaussian 
distribution with a standard deviation given by 

where u is the laser frequency, 7(^0) — w^2I p /F(t Q ) the Keldysh parameter [8] at 
ionization time to, F(t ) is the field strength at ionization time and I p is the ionization 
potential of the atom. 

In [9] experimental evidence was presented for the existence of the momentum 
distribution parallel to the electric field at the tunnel exit. We follow up on this result 
by presenting a detailed analysis of semiclassical simulations and their comparison to 
experimental results using a newly developed method of elliptical integration, which is 
robust at all ellipticities e. Other methods, such as radial integration, work well at high 
ellipticities, but break down for low e, mixing transverse and longitudinal components of 
the momentum spread. On the other hand, projection onto the main axis of polarization 
works well at low e, but breaks down at higher e. 

The elliptical integration method confirms in agreement with [9], that an initial 
longitudinal momentum spread at the tunnel exit fits in well with our experimental 
results. This initial spread was found to be larger than the transverse spread given by 




2. Semiclassical model 

After the quantum mechanical tunnel ionization, the interaction of the freed electron 
with the laser field is considered classically. The electric field of the pulse can be written 

as 

F 

F (t) = (cos(o;t + (f C Eo)x + e sm(ut + y? C Eo)z/) /(*), (2) 

Vl + e 

where f(t) is the pulse envelope, the x-axis is the major axis of the polarization ellipse 
and the y-axis is the minor axis of the polarization ellipse. The momentum spread of the 
electron wave packet at the tunnel exit transverse to the electric field follows from the 
ionization rate calculation. Formally, the kinetic energy from the transverse momentum 
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adds to the ionization potential [10]. In a quasistatic picture, the ionization probability 
at the tunnel exit can thus be expressed as follows [TO] . 



P(, ( . 0. , v ,,) * exp ( - 2( ^ p)3/2 ) exp (- V ^S) . (3) 



3F / V 2aj 

where the longitudinal momentum at the tunnel exit parallel to the electric field, v x , is 
typically assumed to be zero. The transverse momentum spread is given by [5J [7] 

2 _ F 

a± ~ 2(2/ p )V2 ~ 2 7 ' 1 } 

and the exponential was expanded to give two separate terms for the ionization potential 
and the transverse momentum distribution. This transverse spread is valid directly at 
the tunnel exit point as well as at the detector after propagation in the laser field, since 
the field exerts virtually no net force in the transverse direction and the influence of the 
ion Coulomb field is neglected after ionization in SEA [10] . 

On the other hand, there is no clear way of calculating the longitudinal momentum 
spread at the exit point from tunnelling models. But, a theoretical estimate for the 
longitudinal momentum spread at the detector [TTj 



final _ / 3CJ 

a W -y 2 7 3 (l-e 2 )' () 

can be calculated, assuming zero initial longitudinal momentum. This spread is due to 
the different phases of the laser field at ionization time. Electrons ionized before or after 
a peak of the field feel a net force due to the field, which leads to a spread of longitudinal 
momentum acquired during propagation in the laser field. This is verified in figure [TJ 
where ^ is compared to the results from a simulation (find details in section [3]) showing 
very good agreement for e not at extreme values. 

The middle optical cycle has the maximum peak field strength, which decreases for 
optical cycles outside the pulse centre. Averaging the variances from (J5]) with weights 
given by the ionization probability at the peak of each optical cycle P(p) from (J9]) yields 
a slightly smaller overall longitudinal spread 

final/'-\\2\ l / 2 



T.p(p)K M) 



3. Semiclassical Simulation 



All simulations use atomic units. A Classical Trajectory Monte Carlo (CTMC) 
simulation was performed to investigate how the standard assumption of crip 1 * 18,1 = for 
the longitudinal momentum has to be adapted to fit the experiment. For varying initial 
momentum spread at the exit point, the final momentum distribution was calculated. 
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Figure 1. Final longitudinal momentum spread. The analytical formula ^ is 
drawn as a black solid line in both plots, (a) Red x show values obtained from Monte 
Carlo simulations where the ionization events were restricted to the optical cycle at 
the peak of the pulse, (b) The black dotted curve is the averaged analytical formula 
Red x show values obtained from Monte Carlo simulations where ionizations are 
spread over the pulse. 



The simulations follows the TIPIS model (Tunnel Ionization in Parabolic 
coordinates with Induced dipole and Stark shift) [3j. The Stark shifted ionization 
potential is given by 

1 



J p = I p (F(t )) 



( P,0 + 2 ( a N 



ai) F(t c 



(7) 



where % and ai denote the polarizability of the atom and the ion respectively. The 
equation of motion, including the induced dipole in the ion due to the laser field 



r{t) = -F{t)-V[V{r,t)\, V{r,t) 



F r 



ar 



(8) 



is solved numerically to compute the trajectory of electrons post ionization. The 
simulation uses a cos 2 laser pulse centred about t — with a complete pulse duration of 
90 fs, a wavelength of 788 nm and Fq = 0.15 a.u.. These parameters are chosen to match 
the experiment in [9]. Ensembles of 3 x 10 5 events are created for each ellipticity and 
longitudinal momentum spread combination. All electrons have a random exit time, 
with probability weighted by jH] 



P(to) oc exp 



(9) 



3F(t ) 

and an instantaneous tunnelling time is assumed. For the initial condition the TIPIS 
model [3J is used, which calculates the tunnel exit point by solving the cubic equation 
that gives the potential in the parabolic coordinate 77, where r\ corresponds to the 
direction of tunnelling. For 77 ^> 1 (corresponding to exit radius r e > 5 a.u., a condition 
satisfied by present day strong field ionization experiments), the cubic potential in the 
direction of tunnelling is well-approximated by a quadratic, since the terms oc l/r? 2 in 
the potential described in [3] can be neglected. This quadratic approximation yields 



Ve 

2 



^ 2 F(t c 



2F{t ) 



(10) 
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for the tunnel exit radius, while 0% is given by 



02 



2L, 



Figure [2] shows that the exit radius given by the full solution of the potential in parabolic 
coordinates and the quadratic approximation show excellent agreement over a wide range 
of electric field strengths. Additionally, the tunnel exit coordinates given by (10) agree 




Figure 2. Tunnel exit radius. The exit radius calculated from the quadratic 
approximation given by ( |10[ ) (solid line) is compared to the full cubic solution in 
parabolic coordinates (dots). 



within 2.2% with non-adiabatic theory [12] for our experimental parameters (see section 
[5]). Figure [3] shows a comparison of the tunnel exit points predicted by the TIPIS model 
with non-adiabatic, gamma-dependent values given in |12j . 




Figure 3. Non-adiabatic tunnel exit radius. The exit radius as calculated by the 
TIPIS model ( 10 1 (red solid line) is compared to the values found using non-adiabatic 
theory [T^] (blue dashed line) over the complete range of ellipticity. 



To approximate the ion potential ([8]), a soft core potential [I] with a 2 = 0.1 a.u. is 
implemented. 

At any instant, the direction of the laser field defines a coordinate system with 
basis 6_L,i p , 6±, op }, parallel to the field, orthogonal to the field but in the plane of 
polarization, and orthogonal to the plane of polarization. For both y™ 1 *^ 1 and v™ 1 ^ , 
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Gaussian distributed values are generated independently using the standard deviation 
given by Q, resulting in the transverse momentum distribution 



P{v± ,ip> u J-,op) oc exp 



.initial 2 
-L.ip ' 



2a 2 



exp 



.initial^ 

-L,op i 



2a 2 



(12) 



Figure [4] shows momentum distributions calculated from the simulation. As 
expected, with increasing ellipticity the two centres of distribution move to higher 
transverse momentum and the distribution approaches a doughnut shape. For the 
calculation of the distribution, Rydberg electrons jH [13] and electrons which came closer 
than 5a.u. during the pulse are discarded, because at these length scales the classical 
model fails in reconstructing quantum mechanical interactions between the electron and 
the ion. 
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Figure 4. Electron momentum distribution scan over ellipticity. For six 

different cllipticities e, the momentum distribution of ionized electrons projected onto 
the plane of polarization is shown, where x is the major axis of polarization and y 
the minor axis. The slight anticlockwise tilt of the centres is due to the Coulomb 
interaction with the parent ion j3] . 



4. Analysis of final longitudinal momentum distribution 

Here we describe three different methods for analysing the final longitudinal momentum 
spread: 1) radial integration, 2) projection onto the major axis of polarization and 3) 
elliptical segment integration. The intent of all three integration methods is to integrate 
over the transverse momenta, so that the resulting one-dimensional distributions are 
functions of the longitudinal spread only. The first two methods break down at low and 
high ellipticities, respectively, while elliptical integration holds over the entire ellipticity 



Longitudinal momentum spread after tunnel ionization 



7 



range |e| < 1. At fully circular polarization e = ±1, the momentum distribution is 
isotropic and no information about the longitudinal spread can be extracted. This is 
mirrored in the fact that the analytical formula ^ diverges as |e| approaches 1. 

4-1. Radial integration 

Previously, the two dimensional momentum distribution in the plane of polarization was 
integrated radially, and the resulting angular distribution was compared [9] . The radial 
integration is easy to implement and effective for high ellipticity, where it focuses on 
the spread parallel to the field, integrating over the transverse spread. With decreasing 
ellipticity, however, the radial integration mixes transverse and longitudinal components 
more and more as the centre of distribution moves closer to the origin (see figure [4]). This 
causes problems with reliability in the fitting process, and the angular distribution starts 
to develop a spurious double peak structure, where there is only one peak in the two- 



dimensional distribution (see for example figure 12). Therefore, the radial integration 



method is only reliable at higher ellipticities, increasing in accuracy as e increases. 
4-2. Projection onto the major axis of polarization 

As an alternative, the momentum distribution can be integrated over both transverse 
momentum distributions (out of plane of polarization z and minor polarization axis y), 
thus creating a momentum distribution projected onto the major axis of polarization, x. 
This method works well for small ellipticities, where the momentum distribution is only 
slightly curved. However, with increasing e, the accuracy of this projection breaks down, 
as it begins to integrate over the longitudinal spread, particularly at higher absolute 
values of v x (see figure lil). 



4-3. Elliptical segment integration 

It is desirable to have a single technique for extracting the longitudinal spread applicable 
over the whole range of ellipticities. It must take into account that the momentum 
distribution is centred on an ellipse, with eccentricity given by e. To analyse the 
longitudinal momentum distribution, an integration over lines perpendicular to this 
ellipse has to be performed. The momentum distribution as a function of angle is then 
given by, 

f y + 

P(4>) = Vl + e 2 tan 2 / f(x - etan0(y - y ), v) dy, (13) 

Jy- 

y ± =y ±^(l + e 2 tan 2 (f ) )- 1 ^ 

= — /t— t — 9 sin 0, y = — t== cos 0, 
ojyl + e 2 ojyl + e 2 

where is counted from the y axis counter clockwise and f(x,y) = f(v x ,v y ) is 
the momentum probability density in the plane of polarization. In the numerical 
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approximation, this corresponds to binning all events into the ellipse segments pictured 
in figure [5] and weighting the results with the inverse area of the corresponding segment. 



e = 0.6, F = 0.151 a.u., w = 0.058 a.u. 
-1 

-2 



-2 2 

v x / a.u. 

Figure 5. Elliptical segments. The final momentum distributions lie around the 
ellipse shown as a dotted line. The plotted segments are orthogonal to the elliptic 
curve, denoting the segments for elliptical binning of the data to study the longitudinal 
distribution parallel to the ellipse. 

The longitudinal momentum spread at the detector can easily be extracted from 
this technique. Fitting a Gaussian to the resulting momentum distribution yields an 
angle of highest probability (fi m and the standard deviation in angle a^. From this, the 
longitudinal momentum spread is given by 

1 F n 



4 nal = - 7 =^= / a/cos 2 <P + e 2 sin 2 <p d<p. (14) 

5. Experimental Setup 

A COLd Target Recoil Ion Momentum Spectroscopy (COLTRIMS) [H] setup measures 
the ion momentum of Helium ions, which is the negative of the electron momentum 
due to momentum conservation. In the COLTRIMS setup, the fragments created in the 
interaction region are guided towards time and position sensitive detectors by constant 
electric and magnetic fields. The raw data, consisting of the time-of-flight and positions 
of impact of ions and electrons, is used to calibrate the electric and magnetic fields. 
Then, the electron and ion momenta are calculated by solving the equations of motion. 

The momentum resolution is 0.1 a.u. in time-of-flight direction and estimated to 
0.9 a.u. in gas jet direction, mainly determined by thermal spread of the gas jet. The 
momentum resolution is better in the lower half of the detector plate, therefore we 
consider only ions that impact in the lower part of the detector plate and mirror the 
distribution onto the upper half. This makes the spectra perfectly point-symmetric. 

A Titanium: Sapphire based laser system produces short intense laser pulses with 
a pulse duration of 33 fs (FWHM) at a central wavelength of 788 nm. The pulse 
field can be approximated by (J5]). The carrier-envelope-offset (CEO) phase y^CEO [15] 
is not stabilized. The peak intensity / = Fq is estimated to be 0.8 x 10 15 W/cm 2 by 



Longitudinal momentum spread after tunnel ionization 



9 



matching a Monte Carlo simulation described in [9] to the data regarding the momentum 
distribution along the y-axis [16] (the atomic unit of intensity is 3.509 x 10 16 W/cm 2 ). 
The Keldysh parameter depends on the ellipticity and ranges between 7 = 0.51 for e = 
and 7 = 0.73 for e = 1. 

The polarimetry of the experiment has been described elsewhere [3]. Using a 
broadband quarter-wave plate, ellipticities up to e ± 0.93 can be achieved. While 
recording the COLTRIMS data, the quarter-wave plate is rotated continuously by a 
motorized rotary stage and the angle is read out and tagged to the dataset for each 
laser pulse. The angular orientation of the polarization ellipse is calculated for each 
measured ion, and in the presentation of the data the x-axis designates always the 
major polarization axis rather than an axis fixed in laboratory space. The calculation 
of the ellipticity allows generating ellipticity-resolved spectra with a high resolution. 

6. Measurements 

Figure [6] shows the momentum distribution projected onto the plane of polarization 
in the case of anticlockwise rotating field. For each ellipticity, recorded events in an 
interval of ±0.025 around the indicated ellipticity were integrated. 

For e > 0, the two main distributions split apart and lengthen with increasing 
ellipticity, to form a near circular distribution at e ~ 1. Results from the clockwise 
rotating field should be mirror symmetric with respect to the y-axis, compared to 
counter-clockwise field. They were recorded as well and included in quantitative 
comparisons to simulation data (see [4]) to reduce systematic errors. 

7. Longitudinal momentum spread 

In order to find the best fitting initial longitudinal momentum spread from the 
simulation, quantitative comparisons between experimental data and simulation results 
were calculated. Both the experimental momentum distribution as well as the resulting 
distributions from the simulation were analysed using all of the above discussed 
techniques. For each ellipticity, a set of simulations with varying initial longitudinal 
momentum spread, step size 0.05 a.u., was calculated. Figure [7] shows an example of a 
comparison between the simulation and the experimental data for the case of ellipticity 
0.5 as calculated using radial integration (a), x projection (b) and elliptical integration 
(c). 
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Figure 6. Ion momentum distribution scan over ellipticity. For six different 
ellipticities e, the momentum distribution of ionized electrons as recorded by the 
COLTRIMS setup is shown, where x is the major axis of polarization and y the 
minor axis. With increasing ellipticity, the two centres of distribution move to higher 
transverse momentum and the distribution approaches a doughnut shape. The slight 
anticlockwise tilt of the centres is due to the Coulomb interaction with the parent ion 
®. 
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Figure 7. Momentum distribution comparison for the top half. Red dots are 
the experimental data for ellipticity e = 0.5, blue solid lines show the results from the 
simulation where initial longitudinal momentum spreads of a.u., 0.4 a.u. and 0.8 
a.u. were used, (a) Comparison using radial integration, (f> counted anticlockwise from 
the x-axis. (b) Comparison using x projection, (c) Comparison using the elliptical 
segments technique, <j) counted anticlockwise from the y-axis. 
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Figure [7] shows how the initial longitudinal spread is extracted by comparing 
semiclassical simulations to experimental data. For example, an initial longitudinal 
momentum spread cri| lltial = 0.4 a.u. results in good agreement between the simulation 
and the experiment. Both lower and higher spreads (exemplified with <T|| lltial = 
Oa.u. and 0.8 a.u.) result in less agreement. The mean square error for each simulation 
with the corresponding experimental data was computed. A quadratic curve was 
fitted through the errors to find the best fitting initial longitudinal momentum spread. 
The minimum point of the fitted curve was accepted as the best initial longitudinal 
momentum spread to reproduce the experiment. Figure [8] shows the calculated 
simulation errors for e = 0.5 with the elliptical integration. The confidence interval 
for the minimum parameter of the quadratic fit to the set of errors was taken as a lower 
limit error bar length for the respective technique. 



e = 0.5 




0||/a.u. 



Figure 8. Simulation error depending on initial longitudinal momentum 
spread. Evolution of the mean square errors (blue x) between simulations and 
experimental data for increasing initial longitudinal momentum spreads in the 
simulations, calculated from elliptical segment integration. The black solid line shows 
the quadratic curve fitted through the errors, and the red • its minimum point. The 
red O indicate the confidence interval for the minimum position. 



7.1. Results 

The longitudinal momentum spread at the tunnel exit, calculated using three different 
methods (radial integration, projection and elliptical integration), are shown in figure [9j 
along with the theoretically calculated longitudinal spread in (|5]), and the experimental 
longitudinal spread recorded at the detector. The error bars show confidence intervals 
for a confidence level of 0.98 for the extraction of the minimum point in the quadratic 
fit to the simulation-experiment errors. 

For values of e between 0.2 and 0.9, the analysis was performed. As mentioned 
before, only when the polarization |e| < 1 breaks the rotation symmetry, the longitudinal 
spread can be studied. For e < 0.2, there are too many electrons which come closer than 
5 a.u. to the parent ion and have to be discarded in the simulation, such that a direct 
comparison of the calculated distribution to the measured data is not reasonable. 
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Figure 9. Comparison of initial and final longitudinal momentum spreads. 

Black O are the longitudinal momentum spread extracted from the experimental data 
using elliptical integration. The black solid line is the analytical formula ([5| for the 
final longitudinal momentum spread. The black □ show the analytical longitudinal 
spread averaged over the decreasing intensity for all optical cycles in the pulse (JsJ) . 
The best fitting values from the simulation are marked in blue ► slightly shifted to the 
right using the radial integration for e > 0.4, green slightly shifted to the left using 
the x projection, and red ♦ using the elliptical segments technique. All of the error 
bars show the confidence interval for 98% of the corresponding fit parameter. The 
values found show a similar trend as the analytical formula ^ indicated by the red 
dotted line, shifted down to fit the red values. The black dashed line at the bottom is 
the analytical formula for the transverse momentum spread. 



The general trend of the new simulations shows a different behaviour to the one 
given in [9]. The values for cj|| lltial lie almost on a curve given by (5| shifted down to an 
appropriate level. 
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7.2. Effective ionization potential 



To test whether the additional longitudinal spread is caused by different exit points, we 
use an effective ionization potential given by [T7] 

/„, initials 2 



( P,cff 



(15) 



Using this effective ionization potential in the calculation of the tunnel exit makes the 
exit point dependant on the initial transverse momentum of the respective electron at 
the tunnel exit, introducing an additional spread in the final longitudinal momenta 
distribution. To study the significance of this effect, the initial longitudinal momentum 
spread was set to Oa.u. for all simulation runs, and the exit point for each electron 
defined using the above equation. It was found however, that including different exit 
radii in the simulations introduced only a small additional spread that does not account 
for the experimentally measured spread found in [9]. 

A comparison of the results with and without effective ionization potential also 
shows a slight shift in the momentum distribution towards higher v x , see for example 



figure 10 for the case of e = 0.3. 
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Figure 10. Effective ionization potential. Calculating the difference (middle) of 



the momentum distribution with the effective ionization potential ( 15 ) (right) and with 
the Stark shifted ionization potential (left) reveals a shift of the momentum distribution 
in positive v x direction for cllipticity e = 0.3. 



But the shift vanishes for e > 0.5, as demonstrated in the evolution of the differences 



in figure 11 



8. Discussion 



8.1. Double peak structure in angular distribution 

When analysing the velocity distribution observed at the detector in [§] using a radial 
integration, a double peak structure was found in the angular distribution for e < 0.35. 



Only for e > 0.4, does the angular distribution approach a Gaussian, see figure |12| This 
double peak structure is a result of radial integration and can be understood as follows. 

In figure [6] for the case of e = 0.2, black rays from the origin of the coordinate 
system are plotted for two example angles. It can be seen that depending on the angle, 
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Figure 11. Effective ionization potential, ellipticity scan. For lower ellipticity, 
the difference in distributions reveals a shift of the momentum distribution in positive 
v x direction. With increasing ellipticity however, the momentum distributions with or 
without effective ionization potential overlap more, until only the increased transverse 
momentum spread is visible. Yellow or red colours indicate higher probability from 
^p.off, blue tones indicate higher probability from I p (F). 
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Figure 12. Angular momentum distribution. The angular distributions arising 
from radial integration (red points) of the top half distributions shown above in figure 
[6] exhibit a double peak structure for small ellipticity e. Only for e large enough can a 
Gaussian fit (black solid line) reproduce the distribution. The angle <j> is counted from 
the ir-axis anticlockwise. 



a different number of pixels with significant count rate are integrated over. Additionally, 
for angles lying close to the x-axis, the integration goes over longitudinal momentum 
as well. This effect is of course stronger for smaller ellipticities, because the centres of 
distribution are then closer to the origin of the coordinate system, and consequently the 
difference in inclusion of longitudinal momentum is more pronounced. This explains why 
the double peak structure was only found for small enough ellipticity in [9]. In other 
words, the angular momentum distribution calculated from radial integration finds a 
spurious double peak which is due to the elliptic geometry of the momentum distribution, 
and it only appears for ellipticity in a range where radial integration strongly mixes the 
transverse and the longitudinal components of the momentum distribution. 

Taking into account this elliptical geometry of the momentum distribution using 
the above defined elliptical segments results in a well-defined Gaussian distribution for 
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Figure 13. Momentum distribution obtained using elliptical integration 
method. Using elliptical segment integration (red points) on the top half momentum 
distributions shown above in figure [6] yields angular distributions which are well- 
reproduced by Gaussian fits (black solid line) for any ellipticity e. The angle <f> is 
counted from the y-axis anticlockwise. 



8.2. Initial longitudinal momentum spread 

The plotted error bars in figure [9] show the confidence interval from the fitted curve 
through the simulation errors depending on the initial longitudinal momentum spread, 
with a confidence level of 98%. For ellipticity too small, the estimated errors from the 
radial integration were sometimes so large that they would have covered more than the 
range of the plot. The x projection seemed to be more stable even for large ellipticities, 
the error bars did grow considerably, but remained below ±0.06 a.u. at maximum. The 
elliptical segment integration however produces accurate results for all ellipticities. For 
this reason, also the final longitudinal momentum spread from the experimental data 
was extracted using this technique. 

The analytical formula ^ does not take into account that the peak field strength 
for optical cycles which are not at the centre of the pulse is weaker than the overall peak 
strength. Averaging over weaker field strengths in both experiment and simulation leads 
to a smaller longitudinal spread acquired during propagation in the laser field. This in 
turn makes the difference between the theoretical prediction in ^ and the experiment 
in [9] even larger. 

All means of data analysis presented here (radial integration, projection and 
elliptical integration) agree that the fit between the experiment and simulation is best 
if an initial longitudinal momentum spread is included. The values of fx^itiai w j 1 j c ] 1 yield 
the best fit are roughly around 0.4 a.u., and therefore considerably bigger than the 
detector resolution. 

Furthermore, the values of initial longitudinal momentum spread show a strikingly 
similar behaviour to the analytical formula (|5]). Fitting that function through the values 
obtained by elliptical integration results in a reasonable agreement. However, to the best 
of our knowledge there is no physical model that would justify this agreement. 

The ail"* 1111 values which have been previously published in [9] are somewhat higher 
than the values found from this simulation. For small ellipticity, this can be attributed 
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to the radial integration technique, which was used exclusively to analyse the results in 
[9]. Additionally, the initial transverse momentum spread was generated as described 
in [18], which assumes a transverse momentum probability 

/ / initial 2 \ 

P(K nitial |)ocexp -^-^ . (16) 



2ai 



The spread resulting from this distribution is narrower than the one from (12), figure 
14 illustrates the difference. The discrepancy is due to the fact that the probability 



distribution of initial transverse momenta 




1 2 3 

|v | / arb. units 



Figure 14. Initial transverse momentum spread. The absolute value of the 



initial transverse momentum used in [§] given by ( 16 ) (dashed line) and in this new 
simulation obtained from (12) using Q (solid line) show different distributions. The 
average transverse momentum in the new simulation is higher. 



distribution must be transformed as an integral over all t>j_,i p , ^_i_,op which fulfil v\ 
K£ a1 ) 2 + (<% a1 ) 2 , yielding 



/ | „. initial "\ 2 
.initial |\ „, „„„ / \ u . 



P(K ltia l) oc exp y-H^j Z™±- (17) 

The above distribution is valid if 2D CTMC simulation is employed, where the two 
dimensions correspond to the coordinate along the major axis of polarization and the 
perpendicular radial coordinate. 



9. Conclusion 



We present an analysis of the longitudinal momentum distribution using a method of 
elliptical integration. Unlike radial integration (projection) methods, which mix the 
longitudinal and transverse components of the momentum spread at low (high) e, this 
method is robust over the complete range of ellipticity. With the elliptical integration 
method, the longitudinal momentum spreads are well-reproduced by a Gaussian over 
the entire ellipticity range, avoiding the double-peak structure that is observed when 
using radial integration at low e. In agreement with we find that including an initial 
longitudinal momentum spread at the tunnel exit accounts for our experimental results. 
This is in contradiction to standard theoretical assumptions. Further theoretical work 
is necessary to find a physical model that accounts for these results. 
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